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We describe novel purification protocols for bicolorable graph states. The protocols scale efficiently 
for large graph states. We introduce a method of analysis that allows us to derive simple recursion 
relations characterizing their behavior as well as analytical expressions for their thresholds and 
fixed-point behavior. We introduce two purification protocols with high threshold. They can, for 
graph degree 4, tolerate 1% (3%) gate error or 20% (30%) local error. 



I. INTRODUCTION 

The known protocols in quantum information process- 
ing require a certain degree of quantum-mechanical en- 
tanglement to achieve an advantage over their classical 
counterparts. Often, this quantum-mechanical "essence" 
is provided in terms of in-advance-prepared quantum 
states. For example, Bell states are used in a well-known 
protocol for quantum cryptography pj, and schemes 
for multiparty cryptographic tasks using Greenberger- 
Horne-Zeilinger (GHZ) states and other Calderbank- 
Shor-Steane (CSS) states have been devised Q]. Further, 
in quantum computation, multiparticle entangled states 
can be used to streamline the execution of gates and sub- 
circuits via gate teleportation 3], and cluster states rep- 
resent a universal resource for quantum computation by 
local measurements Q. 

In most realistic scenarios the quality of entangled re- 
source states is degraded by the effects of decoherence 
and methods of error detection or correction are required 
to counteract this process. One such method is state pu- 
rification where a (close to) perfect copy of a quantum 
state is distilled out of many imperfect ones. Purifica- 
tion was first described for Bell states 0,IE0 an d subse- 
quentlygeneralized to bicolorable graph states and CSS 
states U S E3 • Recently, a protocol for the purifica- 
tion of W states was presented in ^lj. State purifica- 
tion is used, for example, to establish a perfect quan- 
tum channel between two parties [f| , to efficiently create 
long-range entanglement via quantum repeaters [T^ | or 
to render certain schemes for topological fault-tolerant 
quantum computation universal |13| . 

Imperfect initial states are not the only sources of er- 
ror for realistic state purification. With the exception 
of certain schemes of topological quantum computation 
such as |l3| , errors in the gates for purification also need 
to be taken into account. 

What can we expect to gain from an imperfect purifica- 
tion procedure? In the process of purification the errors 
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of the initial state are replaced by the errors of the puri- 
fying gates. Thus, the amount of error may be reduced 
if the quality of the initial states is low compared to the 
quality of the gates for purification (but above threshold). 
Further, purification can be used to condition the error 
of a quantum state. For example, imperfect Bell-state 
purification can be used to establish a perfectly private 
if imperfect quantum channel |l4j . In a multiparty sce- 
nario, for some protocols the purification gates act locally 
on each copy of the state to purify, resulting in a local or 
close to local error model for the final state. This feature 
attains relevance in the context of fault-tolerant quantum 
computation. Threshold theorems have been established 
for increasingly general types of error including coherent 
and long-range errors jl5l Il6| but there are realistic sce- 
narios in which standard error correction appears to fail 
|l7j . In such a situation, state purification may be used 
to turn the error model into a more benign one. 

The focus of this paper is purification of bicolorable 
graph states by imperfect means, a subject that has pre- 
viously been studied in [jj, [Ha, ■ We are interested in 
the interplay between threshold and overhead. Specif- 
ically, we seek protocols that, (I) work with erroneous 
purification gates, (II) have a high threshold and good 
quality of the output state, (III) scale efficiently, and (IV) 
are analytically tractable. 

Hashing 0, 0, ^| protocols have a high threshold in 
the error of the initial state and require only a mini- 
mal resource overhead, but they break down as soon as 
the purification gates become slightly imperfect [2(j . Re- 
cursive protocols such as [g| also have a high threshold 
for error in the initial states and furthermore work with 
imperfect purification gates, but they are exponentially 
inefficient in the number of particles. 

Our protocols are resistant to initial as well as pu- 
rification errors and are computationally efficient. As 
a bonus, our protocols are analytically tractable for a 
wide class of errors. Specifically, our base protocol de- 
scribed in Sec. IIIII can be analyzed for arbitrary input 
states and general probabilistic Pauli errors in the purifi- 
cation gates. This fact arises through a special locality 
property. So far, the exponential increase of parame- 
ters in the description of n-particle mixed states — even 
mixed stabilizer states — has been found to be an obstacle 
to analytic discussion, and only severely restricted error 
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FIG. 1: Action of MCNOT in the graph basis. The arrows 
represent the direction of syndrome (or Z error) flow (i.e. the 
action of the MCNOT on the stabilizer) 



eigenvalues of all Kj with j e A. Suppose we get -1 at 
the kth qubit. Then we know that either p*; or Uk was 1, 
but we do not have enough information to decide which 
one, so we throw away the states and start again. We 
keep doing this until all measurements are clear. By this 
procedure we correct, to lowest order, errors in the qubits 
of color A. In the next round we interchange the roles 
of colors A and B and so purify the B qubits. We can 
concatenate this procedure to achieve desired levels of pu- 
rity. Because we are post-selecting states on the basis of a 
global measurement outcome, this protocol is inefficient 
for large states. This inefficiency can be addressed by 
using error correction instead of post-selection, to which 
we now turn. 



models have been treated in the literature. 

This paper is organized as follows: in Sec.[n]we briefly 
review the protocol 8] for purification of bicolorable 
graph states. In Sees. IIIII and IIVCI we describe our 
purification protocols and characterize them in terms of 
purification threshold, output quality, and overhead. We 
conclude with a discussion of our results in Section Ivl 

II. BRIEF REVIEW 

Consider a graph G(V,E) with vertex set V and edge 
set E. G(V, E) is bicolorable if V can be partitioned into 
two disjoint subsets A and B such that every edge in E 
connects a vertex in A with a vertex in B. E defines a 
neighborhood relation on elements of V; Af(j) :— {i £ V : 
(i j) E E}. Define the correlation operators 

K j :=X j Zi (1) 

where X, Y, and Z are the Pauli matrices. A graph state 
is a |F|-qubit state \fi) (fi e {0,1}'^') that satisfies the 
eigenvalue equations 

^>) = (-i)^>>,vj = i,...,in (2) 

The states form a basis of the Hilbert space of \V\- 

qubit states called the graph basis. 

We now briefly discuss the post-selection protocol of 
. The protocol works by taking two identical copies 
of a bicolorable graph state and performing multiple 
CNOTs (MCNOT) between them, in a definite pattern 
as illustrated in Fig. ^ Relabeling states in the graph 
basis to reflect the partition into colors A and B (i.e., 
\fjt) = \ti A ,fi B )), the effect of the MCNOT is @ 

\Ha,Vb) \va,vb) i-> \i>-AiPb+vb) Wa + Ha,vb) , 

(3) 

where -I- is elementwise addition modulo 2. Notice that 
information about (jla has been copied into state 2 and 
information about i/b has been copied into state 1. We 
then measure the local observables X and Z on copy 
2, and reconstruct from the measurement outcomes the 



III. THREE-COPY PROTOCOL 

The simplest way to get enough information to perform 
error correction is to do the MCNOT on three copies 
instead of two. The three-copy protocol consists of two 
subprotocols. We use three identical copies of the state 
in each subprotocol. The output of the first subprotocol 
is used as input for the next. Thus, we need nine copies 
to run a single round. Let the three identical copies be 
p (0 \ pW, and p( 2 ). Subprotocol 1 (PI): 

i. Partition the graph into two colors A and B (V = 
V A U V b and V A D V B = 0). 

ii. Perform the MCNOT between copies and 
and and such that information about 
qubits of color A flows from p(°) — > and 
p(o) _> p(2) As a side effect information about 
B will flow from pW, p( 2 ) -> p(°). See Fig. 
below. 

iii. Measure qubits of color A in the X basis and qubits 
of color B in the Z basis in states p^ and . This 
is a measurement of Kj for j £ A. If the measure- 
ment of Kj gives +1 (—1) we get a syndrome of 
(1). Thus, for each j e A we have two bits of 

j (1) , (2) 

syndrome and er- . 

cr (1) -a (2) 

iv. Apply the correction YijeA^j 3 3 to p(°). 

For subprotocol P2 the roles colors A and B are inter- 
changed. 

First, we will analyze this protocol with ideal CNOT 
gates. This will allow us to derive simple closed-form re- 
cursion relations characterizing the behavior of the pro- 
tocol, as well as analytical estimates of the threshold and 
efficiency. In Sec. lTlIBI we generalize to noisv gates. The 
analysis is restricted to density matrices that are diagonal 
in the graph basis (i.e., probabilistic mixtures of graph 
states). At the end of Sec. IIII Bl we will show that our 
results are valid for arbitrary density matrices. 
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A. Ideal gates 



Equation |21 implies that the effect of the MCNOT on 
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Equation [21 implies that the effect of the correction is 
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where o~j := cr^ • "j 2,1 . By measuring p^ and p'- 2 -', we 
get two bits of syndrome for each qubit of color A in p^> . 
The syndrome is conclusive; it allows us to identify, to 
lowest order in the error probability on which state the 
error occurred. We can thus do error correction instead 
of post-selection. This will make the protocol scale effi- 
ciently in the size of the states. The price is a reduction 
of the threshold value. 

We now derive a recursion relation for the expectation 
values (Kj), j s 1, . . . , N. They yield a necessary and 
sufficient condition for purification. For the moment we 
assume that the initial state p is diagonal in the graph 
basis-i.e., p is a probabilistic mixture. It is then safe 
to consider error probabilities. This assumption is not 
necessary, however. It is removed in Sec. IIII Bl Define 
Pj{p) as the probability to find the eigenvalue —I in the 
measurement of Kj on p as 



J2) 



Pj(p) :=Tr 



I - K* 



I - (Kj) 



(7) 



Consider subprotocol PI. In order to analyze this pro- 
tocol we make use of the fact that the error correction 
operation is local. It only uses information about (Kj) in 
each copy to apply a correction to the j'th qubit in p^ . 
Thus, (Kj) should have nice decoupled recursion rela- 
tions. We will later derive the recursion relations for the 
expectation value of arbitrary stabilizer elements, which 
in general are more complex. 



First consider qubits of color B. From Eq. (jJJ, p 



(0) 



(0) 



,(2) 



Since our copies are identical, we have 

PAp m ) = Pj(.P {1) ) = Pj(P {2) ) = Pj- Then, P i * Pf + 
3P,(1 



P 3 ?- 



In terms of expectation values, 
(Kj)' = (Kj)\ 



(8) 



Under concatenation of PI with itself, qubits of color B 
are polluted with (Kj) / ) — * (Kj) j = 0. 

Turning our attention to qubits of color A we note that 
error correction fails if pj = 1 for more than one copy. 



Thus, Pj 
values 



P? 



3P?(1 — Pj). In terms of expectation 
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3-(^) 2 {Kj) 



(9) 
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FIG. 2: Recurrence curves for the three-copy protocol. These 
simple curves fully encapsulate the behavior of the protocol 
with ideal gates. The point of intersection with (Kj) = (Kj) 
gives the threshold. If the gates are too noisy, the protocol 
breaks down, as indicated by the lowest curve. 



Under concatenation of PI with itself, qubits of color A 



are purified with (Kj) m — > (Kj)\o\/o\ = 1 - 

Subprotocol P2 is identical to PI except that the roles 
of A and B are interchanged and the three copies are 
the output states from running Pf 3 times. The three- 
copy protocol is the composition of P2 with Pf. Let 
P = P2oPf; then Eqs. © and JjJ) imply that under the 
action of P 



(Kj)' 




(10) 



The recursion relations If 0|) have, for each color, a 
unique repulsive fixed point in the interval (0, f ) which 
separates the basins of attraction for the trivial fixed 
point at and the nontrivial fixed point at 1 (See Fig.|5J|. 
The upper fixed point corresponds to the perfect graph 
state. Thus, the stated protocol purifies a graph state if 
and only if 



(Kj) > 0.7297 for all j in A 
(Kj) > 0.9003 for all j in B. 



(11) 



We can compare these thresholds to the thresholds for 
the post-selection protocol of 8]. For this protocol, it is 
not known how to derive a threshold for general noise 
or even probabilistic Pauli noise. However, for the par- 
ticular case where only independent local phase flip er- 
rors are assumed for the initial states, recursion relations 
can be derived even for post-selection. Then, the PI 
(post-selection) recursion relation for (Kj) with j S B is 



4 



(K 3 )' = (Kj) 2 and for j g A is (X,)' - T f§^. The 

resulting threshold values are (Kj) th — 0.2956 for j £ A 
and (Kj) th = 0.5437 for j g B. 

Returning to our protocol, it is possible to derive re- 
cursion relations for the expectation values of arbitrary 
stabilizer elements. They are not in general decoupled, 
but there is still a notion of locality. The generalized 
relation allows us to compute the recursion relations for 
stabilizers with small support efficiently Define 

Va \Vb\ 

K a , b .= l[K^l[K b /, (12) 

i=i j=i 

I 



where a g {0, 1} |Va| and b £ {0,1} |1/b| . The factors in 
the first product are the stabilizer generators for qubits of 
color A, while those in the second product are for qubits 
of color B. Then (see Appendix under the action of 
subprotocol PI, 



(Ka,b) 



1 

2H 



E 

oi ,a 2 <?Ca 



(13) 



where / < g iff / 3 = whenever gj — 0. Equations 
and |JS} are special cases for (K ai b) = (Kj) with j £ A,B 
respectively. An interesting feature of this equation is 
that it relates a correlator of weight w = \a\ + \b\ to 
correlators of weight no more than w. This makes it 
feasible to calculate the recursion relations for correlators 
of small weight. 

In order to discuss the behavior of this protocol under 
concatenation with itself, it is useful to switch back to 
probability variables. Then Eq. 1|10[) implies that if the 
protocol is concatenated with itself k times, 



Pj(p(o))V 



P 



th 



(14) 



where P t h is the threshold error probability. The k- 
concatenated protocol requires 3 2fe identical copies, thus, 
the protocol is exponentially efficient under concatena- 
tion. The reduction of error, Eq. Ijl4(l , is not conditioned 
on a particular post-selected syndrome. The overhead in 
number of required initial states is independent of the 
size N of the graph state. We conclude that under con- 
catenation the protocol reaches the reference state |0) 
with efficient use of resources. Contrarily, for the post- 
selection protocol the overhead acquires a dependence 
exp(aTV), with some a > 0, due to post-selection of a 
particular syndrome. 



B. Noisy gates 

Now we investigate what happens to this protocol 
when the CNOT gates themselves are noisy. In the three- 
copy protocol CNOT gates act on the same qubit in two 



states p( m ) and 



We model a noisy two-qubit gate 



channel [i.e., the ST/(4)-invariant channel] 



T« :=(l-pa)[i] + g X; [D^ m) ®D 



(fe, 



(15) 



where D i: j £ {I, X, Y, Z} and k is the qubit index. 
p( fc .™) acts on the fcth qubit of p( m \ The Z gates ap- 
plied in the error-correction steps and the measurement 
of the syndrome are assumed to be noiseless. This is 
natural since the Pauli phase flips Z may be omitted as 
physical operations and instead accounted for in the clas- 
sical syndrome processing. We will include the effect of 
measurement errors in the analysis when we consider the 
more sophisticated protocols, which have higher thresh- 
olds than the three-copy protocol. If we consider the ef- 
fect of T (fc) only on (Kj) in state p(°\ then using Eq. @ 
we can reduce the noise to an effective error. For every 
keV :k£ Af(j) U {j} 



1 



^2 

2 



[1] + ^ 



(16) 



as an ideal gate followed by the two-qubit depolarizing 



If fe ^ Af(j) U {j}, then T^ j) is just the identity map. 
Since every error channel commutes with every CNOT, 
we can model the noisy MCNOT as the ideal MCNOT 
followed by |V| noise channels. 

The error channel Eq. I|15f) is local (i.e., it acts only 
on qubit fe in p( m ) and p^- n '). Also the error operators 
are Pauli operators, which map graph states to graph 
states, keeping p diagonal in the graph basis. Thus we 
can expect the noisy recursion relations to have the same 
form as Eq. I)l()|l. Considering only subprotocol PI, the 
jth qubit in p^ is affected by 2(d+l) error channels. For 
simplicity, we assume all vertices of the graph have the 
same degree d. If this is not the case, then there would be 
a different set of recursion relations for each degree. We 
can then choose d to be the maximum degree, in which 
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FIG. 3: The MCNOT for subprotocol Pi in (a) The 3-copy protocol, (b) The band-aid protocol and (c) The conditional bandaid 
protocol. The dotted lines in c) indicate that the bandaids are applied only if there is an ambiguous syndrome at that location. 
Here we show graphs of degree 2, but these protocols can be applied to graphs of any degree. 



case the recursion relations will be lower bounds for all 
other degrees. The total probability that the jth qubit 

is flipped by an error is 1 ^ p ^ . Thus, for qubits 

of color B, 

{Kj}' = o? (Kjf , (17) 

where a = (1 — P2)^ +1 ^- 

The situation is a little more complex for qubits of color 
A as the error in the MCNOT between p<°) and pW is 
propagated by the MCNOT between p^ and p^ (see 
Fig. OUa)). However, the form of the recursion relation 
remains the same. We get 

(Ki)' = Y (2 + a- 1 -(K j f){K j ). (18) 

For a derivation see Appendix IA 21 Composing subpro- 
tocols PI and P2 we get the recursion relations for the 
three-copy protocol with noisy gates 



(Kj)' = 




(19) 

Here, qubits of color A behave worse. Solving the re- 
cursion relations for fixed points, we find that there are 
two non-trivial positive fixed points (see Appendix^ for 
a > 0.9902. Consider the interval [0, 1]. It has at most 
three fixed points = /o < /i < fa < 1. /o and are 
attractive while fi is repulsive. Thus will be a stable 
fixed point for a > 0.9902 and {Kj) initM > f\. This gives 
a threshold for the noise affecting the gates that scales 
inversely proportional to the graph degree d, 



Pth 



9.8 x 10- 
d+1 



Specifically for degrees 2 and 4 we obtain 



Pth 



0.328 % for d = 2, 
0.197 % for d = 4. 



(20) 



(21) 



This is a rather low value, but it will be substantially 
improved when we consider more sophisticated protocols. 

We now show that the recursion relations Eq. (|19fl are 
valid regardless of whether or not the considered states 
are diagonal in the graph basis. To see this, let us define a 
depolarization operator T> which converts an arbitrary n- 
qubit mixed state p into an n-qubit mixed state pd = Dp 
that is diagonal in the graph basis. V takes the form 



V 



n 



[I] + [Ka, ] 




1} + [Ko,b] 



(22) 



where a and b are vectors in a basis of {0, l}!^- 4 and 
{0, l}l yA respectively. 

We only consider PI, the first round of the protocol. 
It is associated with a transformation PI : p — ► p' = 
R (/9® 3 ). R and V commute-i.e., 

R((Vpf 3 ) =Poi?(p® 3 ), (23) 

for any p. For a proof see Appendix El 
Consider a recursion relation of the form 

(K a Mp' D ))=f a . b ({(Ki.j(PD))}), (24) 

with fa t, some function depending on a, b as in Eq. I|13|l . 

Now, 



(K a , b (p' D )) = Tr 



K a , b R ({Vpf 
= Tr [K a<b V o R ^ 
= Tr [V^K atb )p'] 

= (Ka, b (p'))- 

Similarly, (Kij(pu)) — (Kij(p)), such that 

(K a , b (p'))=fa,b{{(K i , j (p))})- 



[by Eq. Qg) ] 

(trace cyclicity) 
(Pt = V) 



(25) 



Thus, a recursion relation of the form of Eq. (|24|) such as 
Eq. I|19|) holds for all states p and not just for diagonal 
states pd = T>p. 
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IV. IMPROVED PROTOCOLS 

A. Error model 

In the following, we consider a scenario where graph 
states are created locally from product states, then dis- 
tributed to several parties and subsequently purified. Er- 
rors occur in each of these steps-specifically, the following 

• There is a two-qubit error T, Eq. I|15jl. associated 
with each controlled-PHASE gate in the creation 
of the graph state, with probability p2- 

• A local depolarizing error with probability p\ oc- 
curs on each graph state qubit during transmission. 

• Every CNOT gate used in purification carries a 
two-qubit error, Eq. (|15fl . with error probability 
P2- Every measurement is modeled by a one-qubit 
depolarizing channel with error probability p 2 fol- 
lowed by a perfect measurement. 



B. Bandaid protocol 

In order to raise the threshold of the three-copy pro- 
tocol, we will try to combine the strategies of error cor- 
rection and post-selection (which has a higher thresh- 
old). One way to do this is to use small highly puri- 
fied GHZ states-i.e., bandaids, to purify the graph one 
vertex at a time. The usual MCNOT is performed be- 
tween the bandaid and the large graph state as shown 
in Fig. OJb) . This copies information about the central 
vertex into the bandaid which is then measured to give 
a syndrome. Since the bandaid is highly purified (for ex- 
ample, by post-selection), it does not pollute the large 
state much. It is important to note that the error correc- 
tion is still local, and we expect the recursion relations 
to be decoupled as in the case of the three-copy protocol. 

The bandaid protocol also has two subprotocols. The 
first one PI is the following. 

i. Partition the graph into two colors A and B (V = 
V A U V b and Va D V b = 0). 

ii. The bandaids are placed over the large state such 
that each central qubit of the bandaid is over a 
vertex of qubit A for all qubits of color A. Perform 
the MCNOT as shown in Fig.^b). 

iii. Measure the central qubit of each bandaid in the X 
basis and the other qubits in the Z basis. For each 
bandaid multiply the measured eigenvalues. If the 
product is (-1) +1 then the syndrome bit Uj is (1) 
0. 

iv. Apply the correction IX/gA Z°j 3 to the large state. 

P2 is the same as PI, with the roles of colors A and B 
reversed. 



Consider subprotocol PI. For qubits of color B the ar- 
gument is very similar to the three-copy protocol, except 
that each qubit is affected by two gates from each of d 
bandaids. Thus, 



(Kj)' = (l-p 2 ) M (K 3 ){K 3 ) 



b • 



(26) 



where (Kj), is the constant initial purity of the bandaid. 

For qubits of color A, first suppose that the CNOT 
gates are ideal. Then, a simple transfer of purity occurs: 



(Kj)' = (K j ) l 



(27) 



If the gates are noisy, Eq. I|27|l is multiplied by a noise 
factor of the form (1 — p2) jt ^ as in the case of the three- 
copy protocol. There is a subtlety involving the temporal 
ordering of the bandaids. The bandaids do not all com- 
mute with each other. There are 1 + d(d — 1) bandaids 
that affect qubit j. One of them is the bandaid that 
is used to purify the qubit. On average k = d ^ d ~^ f 
the rest will be applied before the purifying one. Any 
effect from the k prior bandaids will be erased by the 
purifying bandaid [see Eq. (|27Jl ]. The purifying bandaid 
has d + 1 noisy CNOTs affecting (Kj)', since the noisy 
MCNOT is modeled as an ideal MCNOT followed by 
noise, no information about the noise is propagated to 
the bandaid. Thus, the noise will commute with the er- 
ror correction procedure. Since a measurement error that 
flips the central qubit of the bandaid will cause us to ap- 
ply the wrong error correction operator, it can also be 
reduced to an effective error as given by Eq. . Thus, 
f{d) = 2{d + 1) + k and we have 



(Kj)' = (i- P2 y 



(Kj) 



Jib ■ 



(28) 



Combining subprotocols PI and P2, we get the recursion 
relations for the bandaid protocol with noisy gates as well 
as noisy measurements 



(Kj 



(I-P2) 2 (Kj) b 

d(d + 3) + 4 

(I-P2) 2 (Kj) b 



for j G A, 
for j E B. 



(29) 



The behavior of qubits of color A is worse, and we will 
use their purity as the final purity of the large state. 

As per our error model in Sec. IIV Al the 
noisy CPHASE, CNOT and measurement gates are 
parametrized by p2- The noisy transmission channel is 
parametrized by p\. For the final result, we need to know 
the quality of the bandaids. We assume that these are 
also created locally, then transmitted and purified. The 
bandaids, however, are of fixed size and may thus be pu- 
rified by the post-selection protocol Q with the higher 
threshold. The output quality of the purified bandaids 
is, to leading order in p 2 , 



(Kj) b = l-(d+l)p 2 , 



(30) 



such that 



l-rf(3d+ll) + 6 
(Kj) = 1 p 2 , (31) 
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FIG. 4: Trade-off curves for the bandaid and conditional 
bandaid protocols (d — 4). The decreasing curves represent 
the breakdown of the post-selection protocol, when there is 
too much error. The increasing curves demarcate the region 
where the final purity of the purified states is higher than the 
purity of the unpurified states. It makes sense to purify in 
the shaded regions. 

for small p 2 [from Eqs. and JSDJ]- As Eq. JJD shows, 
with increasing graph degree the effect of errors in the pu- 
rification process is strongly enhanced. One may there- 
fore ask the question whether it is useful to purify at all 
or whether the transmitted state should be used right 
away. To decide this we compute (Kj) after graph-state 
creation and transmission, 

d(d+l) , , , 

(K j ) = (l-p 2 )^(l-p 1 ) d +\ (32) 

See Appendix [D] for a derivation. We compare this ex- 
pression with Eq. (|29|l and find that there is indeed a 
parameter region where it makes sense to purify. This 
region is displayed for graphs of degree d — 4 in Fig. 01 
It is bounded from above and right by the curve which 
indicates the breakdown of the bandaid purification ac- 
cording to the post-selection protocol [?|. If we use 
post-selection to obtain bandaids of high purity, then 
the threshold of the bandaid protocol for degree d graph 
states equals the threshold for purification of a d+ 1-qubit 
GHZ state with the post-selection protocol [8]. How- 
ever, the output purity of the bandaid protocol is smaller. 
Only above the ascending curve is it advantageous to pu- 
rify. 

C. Conditional bandaid protocol 

In order to correct the cP dependence of the fixed 
point in the bandaid protocol, we will combine it with 



the three-copy protocol. The hybrid protocol, called the 
conditional bandaid protocol, sacrifices in threshold to 
improve the fixed point. The fixed-point behavior, at 
least to linear order in gate noise, is almost as good as 
that of the post-selection protocol. 

This protocol proceeds in the same fashion as the 
three-copy protocol, except that two copies are used per 
round, and wherever a measurement of Kj yields eigen- 
value — 1 (i.e., an error), a post-selected bandaid is ap- 
plied to purify qubit j (see Fig.[3J). For small gate noise, 
we expect to have to apply only a few bandaids per round, 
nonetheless, the threshold is set by the qubits to which we 
have to apply bandaids. Locations where a measurement 
of Kj yields 1 are error free to lowest order. Once again, 
we have two subprotocols, Pi and P 2 each purifying a 
different color. 

The analysis is similar to that used in arriving at 
Eq. I|19fl for the three-copy protocol. However, the sit- 
uation is complicated by the fact that the bandaids are 
applied conditioned on the results of measuring . As 
a result, the recursion relations for the one point cor- 
relators are no longer completely decoupled. We can, 
however, find a simple lower bound on them. 

Define (Kj,) to be the minimum purity of the post- 
selected bandaid. It is a constant. For simplicity we 
assume that all qubits in the bandaid have this purity. 
As before, we assume that the graph of the large state is 
translationally invariant, i.e., all vertices have the same 
degree. The definition (3 = (1 — p 2 ) 2 (Kf,) will be useful. 
Consider qubits of color A in subprotocol Pi, then, by a 
derivation similar to Eq. H18|) . 

(Kj)' = | {2a (Kj) + (K b ) - a (K b ) (Kj) 2 ) , (33) 

where a = (1 — p 2 ) d+1 as before. 

So far, we have been exact. Now consider subprotocol 
P 2 . Again focus on qubits of color A. Break P 2 down 
into two steps. In step one, we apply the MCNOT to p^ 
and p^\ It can be readily verified that (Kj) i— > a (Kj) 2 . 
In step two, bandaids conditioned on the measurement 
outcome are applied to qubits of color B. Let y S {0, l} d 
be the measurement results for the neighbors of qubit 
j. A measurement result of one means a bandaid must 
be applied at that location. If a bandaid is applied to 
a neighbor of j, (Kj) is affected by the errors on the 
bandaid, characterized by (K b ) and by two noisy CNOTs. 
Thus (Kj) i — ► 13^ (Kj). Summing over measurement 
outcomes and including step one, we get 

<*i>' = «(E( E^l )(K 3 )\ (34) 

\k=0 \\y\=k J J 

where q y is the probability of measurement outcome y. 
Unfortunately, q y is a function of the general stabilizer 
expectation values (K a ^b), so we will resort to finding a 
lower bound. Since qo = I — S^^o 1v-> we can rewrite the 
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above equation as 

{K^' = a a ( (1 - 9v) + E E ^ I ^ 

\ k>o\ y \=k J 

>a[l-(l-^)E?»] (Kjf, 

using [3 < 1 to arrive at the inequality. 

Now, qo is just the probability that no error is de- 
tected on any of the neighbors of j. Let pj be the 
probability of detecting an error on site j. Then, by 
definition, (Ki) = 1 — 2 X^|y 4 =i %• This implies that 

J2 y ^o1y - J2ieAf(j) ^~~2^ ■ Putting this into the above 
inequality, 

(Kj)' > a (l - |(1 - /3 d ) (1 - (if,))) <i^} 2 , (35) 

where is the purity of qubits of color B from the 
previous round. 

Solving for the fixed point, we get, to leading order in 
gate noise p2, 

(Kj) = 1 - 2(d + l)p2- (36) 

Comparing this to Eq. H3U|) . we see that the fixed- 
point scaling with degree is almost as good as in the 
post-selection protocol. We now apply the conditional 
bandaid protocol to the same situation-of a graph state 
being shared among widely separated parties, as for the 
bandaid protocol. The results for a degree four state 
are plotted in Fig. 0] We see that the threshold (upper) 
curve is worse, whereas the fixed-point (lower) curve is 
better for this protocol, as compared to the bandaid pro- 
tocol. The total purifiable area is smaller, indicating that 
it breaks down faster. In some sense, we have traded 
threshold for fixed point. These conclusions hold for ar- 
bitrary degree, and the curves are independent of the size 
of the state, making this protocol eminently suitable for 
the purification of large bi-colorable graph states. 

V. CONCLUSION AND OUTLOOK 

We have described novel purification protocols for bi- 
colorable graph states and discussed their performance. 
The criteria for our protocols are that they do not break 
down in the presence of small amounts of noise in the 
purification process, that they have a high purification 
threshold and good output quality, scale efficiently, and 
be analytically tractable. 

Our final protocol can, for relevant graph states of de- 
gree 4, tolerate 1% gate or 20% local transmission error. 
These are about 1/3 and 2/3 of the respective values for 
the post-selection protocol UH. However, in contrast 
to this reference protocol, our protocol scales efficiently 
with the graph size. 



All our protocols can be treated analytically. In par- 
ticular, for the three-copy protocol we derive closed, ex- 
act one-dimensional recursion relations in the appropri- 
ate observables, irrespective of the size of the state. 

We would like to comment on the influence of the graph 
degree for the purification threshold. First note that for 
the three-copy protocol of Section III, in the case of per- 
fect purification gates, the recursion relations (|L(J|) are 
completely independent of the graph structure, and so 
are the thresholds . This behavior changes if noise is 
included in the purification. The critical noise level per 
purification gate — at which the protocol breaks down — 
scales inversely proportional with the graph degree. The 
unfavorable dependence on the graph degree is present 
in all three protocols we discuss. Thus, the lesson we 
learn for the case of noisy purification is to beware of 
large graph degrees. Large graph degrees occur, for ex- 
ample, in graphs states corresponding to codewords of 
concatenated CSS codes. 

We would also like to comment on the structure of 
the nontrivial fixed point in our protocols. In the case 
of erroneous purification gates, the nontrivial fixed point 
is not completely specified by the lowest order expecta- 
tion values (Kj) and it remains to be discussed which 
error correlations are removed by the purification pro- 
tocol. As a first result in this direction, for the three- 
copy protocol discussed in Section III we have shown (in 
Appendix I A 3|) that correlations of stabilizer expectation 
values located on non-overlapping supports are not in- 
troduced by the purification procedure if they are absent 
initially. This implies that such correlations are absent in 
all purified states which end up at the same fixed point 
as the perfect state. We show in Appendix [B] that the 
fixed point for two-generator correlations with distinct 
support is unique, which is enough to establish the re- 
sult that all states at the fixed point obey the relation 
(KiKj) = (Ki) (Kj) for such correlations. 

A question of further interest is whether the nontrivial 
fixed point of the protocol is unique at all levels of cor- 
relations. This would imply (Ki + j) = (Ki) (Kj) for all 
correlations with distinct supports. 

Another question of further interest is whether the de- 
scribed or related protocols may be used to boost the 
threshold value for fault-tolerant quantum computation 
HH0 |H i|| based on graph states. 
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coefficient of [I] contributes, which is (— 1) a j 
eluding the stabilier operator, we are left with 

^b(°) bWW) b< 2 ) X 



In- 



^oCi+oW+oW.K 1 )! ( A4 ) 



APPENDIX A: GENERALIZED RECURSION 
RELATIONS 

We now derive the generalized recursion relations 
[Eq. 1(13(1 ] for the three-copy protocol. While the method 
used for this derivation is less intuitive, it yields recursion 
relations for arbitrary stabilizer elements and can handle 
noisy gates easily. 



1. Noiseless gates 

In order to derive Eq. 1)13(1 we work in the stabilizer 
basis. Because p^ is diagonal and the set {(K aib )} 
where a 6 {0, l}' 17 - 4 ', 6 € {0, l}' 17 - 4 forms a complete 
set of observables, we can write an expansion — 

2 \v A \ + \v B \ J2a,b ( K aM) K aib . 

Consider subprotocol Pi, which purifies the A sub- 
graph. The initial state is p^ ®pW ® p^ 2 \ which can be 
rewritten as a sum over a, b of terms of the form 

(^a(°),b<°)) (KaW.bw) ( K a( 2 >M 2) } X 

K a (.0)^0)K a (l) tb (,l)K a (2) tb (2). (Al) 

The protocol is linear, so we track the evolution of each 
term seperately. Performing step (ii) , this term becomes 

(-^oW,6(°)) (-^aW.&m) ( K a( 2 \b' 2 l) x 

K a {0) +a (l) +a {2) b {0)K a[ l) bl a) +b (l)K al 2) b {0) +b (2) . (A2) 

Now consider step (iii). Suppose we get measurement 
outcomes , for the stabilizers in subgraph A on 
copies pd>,p( 2 '. Then the resultant state is given by 
applying the projector 



. \Va\ 

pfcr II + (-D^fXw + (-i)^^ 

(A3) 

All the single-site operators involved commute, so this 
term is a product of stabilizers in b and terms of the 
form 



[I] + (-l)4Vf)(^f) af) 



Here k = 1,2. Discarding p^',p( 2 ', we perform a par- 
tial trace over these systems (recalling that K a b are all 
traceless except Ko,o = I)- In the above term, only the 



where S p ^ q is the Kronecker delta on each component of 
p, q. Note that we must have b^ ' — b^ — b^ 1 or the 
term is zero. 

Now examine the action of the Pauli [Z] operator in 
this basis. [Z]K a , b = ZK a , b Z = — l fe if a ,b, where k = 
iff Z and K a>b commute. Effectively, Z is a diagonal 
matrix with entries ±1. Identical reasoning applies to X 
and Y. This will make it very easy to add gate noise into 
the analysis. It also allows us to say that the net effect of 
the error-correction step iv is to multiply Eq. (|A4|) by a 

factor Of (_l)(A (1 'xA< 2 )).(a'°)+a( 1 '+a< 2 )) ; where {p x ^ = 

Pj ■ qj. To simplify the notation, change the basis to 
a = a(°) + a^ 1 ) + a^ 2 ) , b = b (0) . Then the term becomes 

(_ 1 ^A( 1 )-a< 1 )+A< 2 )-a< 2 ' + (A (1 )xA( 2 ))-a 

S b,bW S bM 2 > 1 2PWI\ Kab ' 

In this notation and ignoring the delta func- 
tions, the original coefficient in Eq. (| A2|> is 

( K a+aW+am,b) ( K aW ,b) ( K aW ,b) ■ We will nOW 

get conditions under which this term contributes to the 
coefficient of K a b . 

Summing over measurement outcomes, the coefficient 
of K ayb is 

( K a+aW + a ( 2 >,b) ( K a W,b) ( K a<- 2 \b) X 

i , -\( 1 ).„.< 1 )+A( 2 '-a( 2 ) + (A< 1 )xA <2 )l 



2 2 \ v a\ 



A(!),A( 2 ) 



The sum can be reexpressed as 



l 



\Va\ 

n ^ 



If a 3 

J 2 ) 



(i) 



= 0, then the jth factor is zero unless a 
0, in which case it is 4. Hence, for the term 



( K a+ a w+ a ( 2 ),b) ( K a w,b) (^a( 2 ) : b) to survive the proce- 
dure, we must have a^\a^ 2 ' <C a. If this holds, then 
an overall factor of 4l v - 4 H a l comes out. If a,j = 1, then 

a straightforward calculation shows that the jth factor 

(i) < 2 ) 

contributes a factor of 2 (—l) a 3 a j . The overall numeri- 
cal factor is thus To get the new value of (K aib ), we 
simply sum over since these and only these will 

contribute to the support of K a b under PI. This gives 
Eq. O- 
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2. Noisy gates 



Behavior of correlations 



Adding noise to the gates requires very little additional 
work. We can rewrite the depolarizing channel on qubit 
j of copy k as 



D 



[P] = ~ ([/] + \z)f ] ) \ ([J] 



P^Px 



If we take two qubits j, k such that Af(j) f)Af(k) — 0, 
then the noise terms on sites in Af(k) U k do not affect 
terms involving j and vice versa. Hence the sum over 
terms in the recursion relation for (Kjk) will factor into 
(Kj){K k ). If initially (KjK k ) = (Kj) (K k ), then the 
three-copy protocol will not generate any new correla- 
tions between these regions. 



It was shown above that [X], [Z] have ±1 on the diagonal. 
Thus writing the noise channel in this form illustrates 
how the noise components act as projectors P%, 



x, 



If a specific ket is affected by noise on site j of copy k, it 
will be an eigenvector of with zero eigenvalue. 

The noise from a CNOT at site j between copies i and 
k is 



E 



(»),(fc) 



(1-P2)+P2(PI'P ( xKPz 



WwpW p( fe ) 



P y x)- (A5) 



If a ket K a b is affected by any of these noise terms (that 
is, if the noise anticommutes with K a ^), it will be pro- 
jected to zero and thus acquire a (1 — pa) multiplier over- 
all. 

The noise from the first MCNOT is E 01 = JJ. £j 0) ' (1) 

and from the second MCNOT is E 02 = HjE^'®. 
Clearly the overall multiplier is independent of the mea- 
surement outcomes, so the analysis for Eq. (|A3I) still 
holds. The recursion relations are then similar in struc- 
ture to Eq. (|13|1 . except that coefficients dependent on 
(1 — P2) are inserted before each term. 

We illustrate this by calculating the recursion relations 
for {Kj). If j £ B, there is no sum in Eq. 1)13(1 . and 
{Kj) — > Ej {Kj) 3 . The only noise terms that anticom- 
mute with Kj [and hence give factors of (1— P2)] are those 
in j U Nj. There are 2(d + 1) of these (since there are 



\2(d+l) 



(Kj) 



two sets of noisy gates), so {Kj) — » (1 — p%Y 
which is Eq. (jTTjl . 

Now suppose j £ A. Let j = (0, . . . , 0, j, 0, . . . , 0). 
Our sum is over a^\a^ £ {0,j}, and b = 0. Since 
we are interested only in (Kj), our effective noise model 
is [X] k 1 — ► [Z]jik £ Nj and [X}j h-> [I]. All other noise 
terms do not affect the state. Then 



(l- P2 )+P2Pr'P. 



(jfi) P 0M) 



(A6) 



A similar replacement holds for Eq 2 - Eqi acts on terms 



E j+aW + a < 2 ) qE a(i) t o and gives a factor of 1 iff j + a 



(i) 



,(2) 



a^' = 0,oW = j = o ffl ,oW = 0, and a factor of 
(1 — P2) d+1 otherwise. 

Performing the MCNOT between and p (2 \ the 



noise channel E$2 acts on the kets K 



j+aW,0 A a( 



K„i2) . which 



gives a factor of 1 iff j + a^ = 0,a^ = and {l-p) d+l 
otherwise. Putting in each of the four cases ,af £ 
{0,1} gives us Eq. flS). 



APPENDIX B: UNIQUENESS OF THE FIXED 
POINT 



Here we show that the three-copy protocol has a unique 
fixed point for stabilizer elements (K a ^) with weight 
w = \a\ + \b\ < 2. The recursion relations for stabi- 
lizer elements of weight w > 1 [see Eq. JT3J] depend 
only on stabilizer elements whose weight is at most w. 
Thus, we can use an inductive argument. If all the stabi- 
lizer elements of weight less than w have reached a fixed 
point, they become constants and then the recursion re- 
lation for elements of weight w will have the same form 
as those for weight one (i.e., they will depend only on 
stabilizer elements of weight w). First consider the case 
when \a\, \b\ < 1. For this case, the three-copy recursion 
relations Eq. (|13|) have the form 

f(z) — az + bz 3 . 
g(z) — cz + dz 3 , 

with a, c > and bd < 0. The presence of noise does 
not change the form of the recursion relations, it only 
multiplies each term by a number between and 1 (see 

Define 



Appendix IA 2|) . Let y — z and x = dy 



p(x) := f(g(z))/z — 1 = &.t 4 — bcx 3 + adx — d. 

The signature of p(x) is 

p(x) :- + +-, 
p(-x) : . 

Then by Descartes' rule of signs |25|. p(x) has at least 
two complex roots. Thus the recursion relation f(g(z)) = 
z has at most two positive fixed points. The recursion 
relation g(f(z)) = z can be analyzed identically. It was 
already argued in Sec. IIIIBl that this means that there is 
a unique attractive fixed point. 

Now consider the case |a| = 2 and |6| =0. The recur- 
sion relations now have the form 

f(z) = az 3 + bz + c 
g(z) = dz 3 . 

It is easily checked that a, c, and d are positive. The 
sign of b is harder to fix, but we note that for there 
to be a fixed point at all, b must be negative. The 
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case f(g(z)) = z is easily analyzed, as above, to show 
that there are at most two positive roots. Let p(z) — 
g(f(z)). To conclude the proof we need two technical 
results, (i) If the smallest support expectation value 
(K a ) has reached its fixed point value (K a ) fp then the 
physically allowed values for (K a+a i) form the interval 
/ = [2(K a ) lp - 1,1]. (ii) f(z) > for all z £ I. 

Proof of (i) (a) z allowed =>• z £ L P = 1- f "' , 

with a ^ a', is a projector, hence (P) > 0. Thus 
z = (K a+a ,) > (K a ) + {K a r) - 1 (*). Evaluate (*) at 
fixed point (K a ) fp . z < 1 is obvious, (b) z £ I => z al- 
lowed: For an initial state of the protocol, interpolate be- 
tween pi = (K a ) {p p ++ + (1 - (K a ) ip )/2 {p+-+p-+) and 
P2 = (K a ) {pP++ + (l-(K a ) (p )p^. (The signs "±" refer 
to the eigenvalues of K a and K a i, respectively.) Proof 
of (ii). Be {K a ) lp , {Kb) {p > and z £ I. Assume as an 
hypothesis f(z) < 0. Apply (*) to the state after appli- 
cation of PI, at the fixed point {K a ) fp , Va £ A. Hence 
> f(z) > 2 (K b ) lp - 1. (Under PI the fixed point value 
(K a ) fp for a £ A is mapped to {Kt) fp for b £ B, assuming 
all vertices have the same degree.) Thus, (Kb) fp < 1/2. 
But then {Kb) lp = 0, which is a contradiction. Hence 
/(*)>0. 

Now, p"(z) = g"(f(z))f'(z) 2 +g'(f(z))f"(z) such that, 
with (ii), p" > for all z £ I. Thus, p(z) is convex on 
/. With (i), I is a single interval such that p(z) and z 
intersect at most twice in /. At most one of these fixed 
points is attractive. 



-. 1 k 




FIG. 5: Creation of a degree (d — 4) bicolorable graph state. 
The figure will have the same local structure for other degrees 
and topologies as long as its edges are d colorable and its 
vertices are bicolorable 

and classical post-processing. K^ b , K^ b commute with 
the Kraus operators in Eq. I|22J) . such that 

Tr AfoD (0) D (1) P (2 ' = Tr V {0) V {1) V {2) oM 

(1,2) (1,2) 

= X>(°) o Tr M. (C4) 

(1,2) 

Eqs. EH, O, and El yield Eq. (gSJ) 

APPENDIX D: CREATION OF A 
BI-COLORABLE GRAPH STATE 



APPENDIX C: THE DEPOLARIZING 
OPERATOR 

In order to prove that the depolarizing operator T> de- 
fined in Eq. (1221) commutes with the evolution opera- 
tor R = Ttm 2) M o £ o U, we note that the protocol 
step PI consists of a unitary part U , an error channel £ 
comprising probabilistic Pauli errors, and a measurement 
Tr( 12 ) M, where M is a projector. U consists of a set of 
transversal CNOT-gates and acts on the stabilizer as 



K (o) 



K 



(2) 



K i0) K (1) K {2) 

a,o O.o' 

K i2) K {0) 

A a, 6^0. fa- 



(CI) 



etc., such that 



Now note that 

w+iK l] m+[Q [i\+[k§1\ 

2 2 2 

U O p(0)p(l)p( 2 ) = p(0)p(l)p(2) Q ^ ( C2 ) 

The operations £>(°)Z)(i):Z?( 2 ) and f commute because 
both are linear combinations of Pauli superoperators, 

£ o V^V^V^ = p(0)p(Dp(2) £> ( C3 ) 

The measurements comprising Tr( 12 ) M are of stabilizer 

operators K^l, K^l on the states p^\ p^ 2 \ respec- 
tively. They are performed via one-qubit measurements 



Here we discuss the noise structure of a bicolorable 
graph that is created using noisy CPHASE gates. The 
noisy gates are modeled as the ideal gate followed by 
two-qubit depolarizing noise as defined in Eq. (|15fl . The 
graph state is created by performing CPHASE gates be- 
tween qubits in the |+) state. The noise structure of 
the final state depends on the temporal ordering of these 
gates. If we assume that the underlying graph has con- 
stant degree d and that its edges are d colorable, then the 
./V-qubit graph state can be created in d time steps with 
Nd CPHASE gates. At each time step all the gates cor- 
responding to edges of a particular color are performed. 
Thus, at every time step t £ {1, ...,d}, each qubit is 
affected by an error channel of the form of Eq. (|15|l . 

We are interested in the value of (Kj), so we focus on 
the neighborhood of qubits j in the larger graph. Since 
the graph is bicolorable, it contains no three cycles and 
one can draw a diagram of the form of Fig.0 The gates 
are represented by both solid as well as dashed lines. 
The noise channels corresponding to the solid lines each 
contribute an effective error T c g as defined in Eq. 1|16|) 
to qubit j. Now consider the qubit k which is a neigh- 
bor of the central qubit j. Each dashed line also con- 
tributes an effective error T e ff to qubit j , but only if the 
CPHASE gate corresponding to the solid line between 
k and j was performed in a previous timestcp. This 
is because errors commute with Kj and Xk errors 
would be propagated by the CPHASE to X^Zj errors, 
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which also commute with Kj. Thus there are a total of 
d ^ 1 " 1 + d = d l d ±ll noise channels affecting the qubit j. 
This gives 

d(d+l) 

(K j ) = (l-p 2 )^. (Dl) 
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